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Abstract A survey was carried out in search for bat 
coronaviruses in an urban maternity roost of about 500 
specimens of two species of insectivorous bats, Molossus 
molossus and Tadarida brasiliensis, in Southern Brazil. 
Twenty-nine out of 150 pooled fecal samples tested posi- 
tive by reverse transcription-PCR contained fragments of 
the RNA-dependent RNA polymerase gene of coronavirus- 
related viruses. The sequences clustered along with bat 
alphacoronaviruses, forming a subcluster within this group. 
Our findings point to the need for risk assessment and 
continued surveillance of coronavirus infections of bats in 
Brazil. 
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Bats (order Chiroptera, suborders Megachiroptera and 
Microchiroptera) are one of the most diverse and widely 
distributed groups of mammals, representing ~20 % of all 
known mammalian species [1]. About a 100 different viruses 
have been identified in bats of different species in Asia, 
Europe, North America and Africa. Therefore, such species 
may be natural reservoirs for a large variety of potentially 
zoonotic RNA viruses, such as lyssaviruses, paramyxovi- 
ruses, Ebola and Marburg viruses as well as the recently 
emerged severe acute respiratory syndrome coronavirus 
(SARS-CoV) [2-6]. A variety of other coronaviruses have 
been detected in many bat species from Asia, including 
specimens of the genus Rhinolophus, which were found to be 
infected with SARS-like CoV. Phylogenetic analyses of 
such viruses revealed that those form a large clade within 
Betacoronavirus genus, along with SARS coronaviruses 
from palm civets and the SARS coronaviruses recovered 
from humans during the 2003 outbreak [7, 8]. These data 
suggested that the agent responsible for the 2002-2003 
pandemic might have originated from bats. In addition, in 
2012, a new human coronavirus (HCoV-EMC), which has 
been associated to clinical disease that resembles SARS, 
emerged in the Middle East. This new virus appears to have 
originated from bats, raising the possibility that HCoV-EMC 
jumped species directly from bats to humans [9]. 

In Brazil, most studies looking for associations between 
bats and viruses have focused on the role for those species as 
reservoirs for rabies virus [10]. However, to date, more than 
160 bat species have been detected in Brazil, comprising 
members of the families Phyllostomidae, Vespertilionidae, 
and Molossidae. It is estimated that at least 40 bat species live 
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in the state of Rio Grande do Sul, Southern Brazil, where the 
predominantly sub-tropical climate seems to favor the 
settlement of such species [11]. 

In view of the potential role that bats may play in the 
transmission of new viral infections to humans and other 
species, this study was set up in search for coronavirus 
genomes in bats from the urban area of Porto Alegre 
(30°01'59"S; 51°13’48”W), a town with about 1.5 million 
inhabitants and capital of the state of Rio Grande do Sul, 
Brazil. With that purpose, coronavirus RNA was searched 
in feces of two species of synanthropic insectivorous bats 
collected in a maternity roost within the urban area of the 
city. 

A maternity roost of bats known to have direct contact 
with people and domestic animals was identified in the 
summer of 2012 in the attic of a residence in the central 
area of Porto Alegre, Southern Brazil. The colony was 
estimated to harbor about 500 bat specimens of insectivo- 
rous bats of two species, velvety free-tailed bats (Molossus 
molossus) and brazilian free-tailed bats (Tadarida brasili- 
ensis). Speciation was confirmed by amplification and 
sequencing of the mitochondrial cytochrome b (cytb) gene 
as described [12]. 

One hundred and fifty fecal samples were collected from 
the attic floor as follows: a plastic film was spread on the 
ground of the attic compartment and fresh droppings were 
collected with clean disposable forks in the following 
night. Each sample consisted of five fecal droppings, which 
were immediately sent to the laboratory and stored at 
—80 °C. The samples were then submitted to total RNA 
extraction with TRIzol (Invitrogen™)., CoV RNA screen- 
ing was performed by reverse transcription-polymerase 
chain reaction (RT-PCR) in a total volume of 25 pL 
reaction using conserved primers for the RNA-dependent 
RNA polymerase gene (forward: 5’-GGTTGGGACTATC 
CTAAGTGTGA-3’ and reverse: 5’-CCATCATCAGATAG 
AATCATCATA-3’). This pair of primers is expected to 
give rise to amplicons of 440 bp [7]. 

The cycling conditions were: 5 min at 94 °C followed by 
35 cycles of 1 min at 94 °C, 1 min at 49 °C and | min at 
72 °C, followed by a final extension time of 5 min at 72 °C. 
Bovine coronavirus (BCoV) RNA was used as a positive 
control to optimize the assay. Standard precautions were 
taken to avoid PCR contamination; blank controls without 
template were included in every set of five RT-PCR assays. 
Five microliters of the PCR products were electrophoresed in 
1.5 % agarose gels and the products visualized on UV light 
after staining with ethidium bromide. The amplicons obtained 
were cloned into pCR® 2.1-TOPO® cloning kit (Invitrogen) 
before being submitted to nucleic acid sequencing. 

Sequencing was performed with the Big Dye Terminator 
Cycle Sequencing Ready Reaction (Applied Biosystems, 
UK) in an ABI-PRISM 3100 Genetic Analyzer (ABI, 
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Foster City, CA), following the manufacturer’s protocol. 
Sequence analyses were performed with the BLAST soft- 
ware [13]. Nucleotide sequences were aligned and com- 
pared to human and animal CoV sequences available at 
GenBank database with the program ClustalX 2.0 [14]. 
Alignments were optimized with the BioEdit Sequence 
Alignment Editor program version 7.0.9 [15]. The protocol 
to generate the phylogenetic trees was selected with the 
program Modeltest 3.7 [16]. Phylogenetic analysis was 
carried out using MEGA 5.0; pairwise genetic distances 
were calculated by the Tamura 3-parameter model and 
phylogenetic trees were constructed using the neighbour- 
joining method. Bootstrap values were determined by 
1,000 replicates to assess the confidence level of each 
branch pattern. 

PCR amplicons with the expected size of the targeted 
region were obtained from 29 out of the 150 (19.33 %) 
pools of bat fecal samples. The nucleotide sequences of 
sixteen randomly selected amplicons were determined and 
submitted to GenBank (accession numbers KC 110770 to 
KC 110785). Genetic analyses provided evidence that the 
viruses circulating in these two species of insectivorous 
bats belong to the genus Alphacoronavirus. When com- 
pared with each other, all the obtained sequences showed a 
high nucleotide and amino acid identity (90.6 to 100 % and 
98 to 100 %, respectively) (supplemental material). The 
RdRp sequences examined here were distantly related 
(<75 % nt identity) to other known alphacoronaviruses. 
The closest bat coronaviruses RdRp sequences found in 
GenBank were the Asian (BtCoV/A633/2005) and North 
American (RM-BtCoV 6 and RM-BatCoV 11) bat coro- 
naviruses (Fig. 1). The percentage of nucleotide similarity 
between the sequences described here and those of Asian 
and North American coronaviruses ranged from 72.4 to 
76 %, whereas at the amino acid level, the similarity ran- 
ged from 74 to 81 % (data not shown). 

During the last two decades, several studies have shown 
that various important human and animal pathogens are of 
bat origin; these species have become targets for several 
surveillance studies aiming the detection of other poten- 
tially pathogenic viruses for humans and other animals. 
The association of these pathogens and possible disease 
outbreaks caused by direct or indirect contact of humans 
with bats stimulated the development of research activities 
on bat-borne viruses. In addition, the advances of molec- 
ular techniques offer opportunities for the discovery of 
novel DNA and RNA bat viruses without the need for virus 
isolation and bat pooled fecal samples being used as source 
for viruses, preventing animal manipulations [17, 18]. 

In our study, we detected RdRp sequences of bat CoV at a 
frequency of 19.33 % in the examined samples; such fre- 
quency is comparable to previous results obtained in similar 
studies from different bat species in other countries (ranging 
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Fig. 1 Phylogenetic 
relationship between the 
nucleotide sequences 
corresponding to a portion of 
440 bp of the RNA-dependent 
RNA polymerase gene 
amplified from insectivorous bat 
species in Porto Alegre, RS, 
Brazil. Sequences from other 
alpha and betacoronaviruses 
were obtained from GenBank 
(accession numbers in the 
figure). The tree was generated 


based on the neighbor joining 
method in the MEGA program. 
The nucleotide sequence of the 
equivalent genome fragment of 
SARS-CoV was included as 
outgroup 
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KC710770 to KC110785| BatCoV/POA/2012 


EF544567.1| Rocky Mountain Bat Coronavirus 3 
DQ648838.1| Bat coronavirus (BtCoV/860/2005) 
GU065403.1| Kenya bat coronavirus BtKY59 
DQ648823.1| Bat coronavirus (BtCoV/A527/2005) 
DQ648825.1| Bat coronavirus (BtCoV/A604/2005) 
EF544563.1| Rocky Mountain Bat Coronavirus 11 
EF544565.1| Rocky Mountain Bat Coronavirus 48 
EF544568.1| Rocky Mountain Bat Coronavirus 6 
54 EF544564.2| Rocky Mountain Bat Coronavirus 27 
JQ765563| Human coronavirus NL63 
53 — YQ023162| Porcine epidemic diarrhea virus 
JQ023161| Porcine epidemic diarrhea virus 

99 | JN825712| Porcine epidemic diarrhea virus 

86 L JQ282909| Porcine epidemic diarrhea virus 
DQ666337.1| Bat coronavirus 1A 
GU065410.1| Kenya bat coronavirus BtKY66 

GU065405.1| Kenya bat coronavirus BtKY61 

DQ648799.1| Bat coronavirus (BtCoV/A1196/2005) 
EU834956.1| Bat coronavirus M.aus/Australia/CoV180/1996 
EF544566.1| Rocky Mountain Bat Coronavirus 65 
GU190237.1| Bat coronavirus 1B BR98-19/BGR/2008 
AY532270.1| Acinonyx jubatus coronavirus 
DQ010921| Feline coronavirus 
AY994055| Feline infectious peritonitis virus 


EU769560| Chinese ferret badger coronavirus 
AF124986| Canine coronavirus 

DQ201447| Transmissible gastroenteritis virus 

DQ443743| Transmissible gastroenteritis virus 

99 'AJ271965| Transmissible gastroenteritis virus 

HM245926| Mink coronavirus strain WD1133 
HM245925| Mink coronavirus strain WD1127 
GU065416.1| Kenya bat coronavirus BtKY72 
DQ648795.1| Bat coronavirus (BtCoV/A1018/2005) 
DQ648803.1| Bat coronavirus (BtCoV/242/2005) 


GU065420.1| Kenya bat coronavirus BtKY76 
GU065398.1| Kenya bat coronavirus BtKY54 
JX174638.1| Bat coronavirus Gabon/292/2009 
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from 2 to 55 %) [19-25]. The findings reported here dem- 
onstrate that the bats in this particular maternity roost are 
infected by coronaviruses, genomes of which may be 
detected in their fecal droppings. Due to methodological 
limitations, we were not able to determine to which bat 
species the bat CoVs-positive samples belonged. Never- 
theless, it becomes clear that such viruses do circulate 
within bats that share this particular roost. Although all 
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AY268070.1| SARS coronavirus Hong Kong 03/2003 


tested samples clustered within Alphacoronavirus genus, 
the RdRp sequences of Brazilian bat CoVs formed a single 
subgroup within it, yet distinct from previously identified 
coronavirus RdRp sequences recovered from a number of 
species. Our samples clustered with another known alpha- 
coronaviruses found in Myotis ricketti BtCoV/A633/2005, 
and Myotis occultus Rocky Mountain bat coronavirus 6 and 
11, reported previously in Asia and North America [18, 24] 
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(Fig. 1). These results show that similar coronaviruses are 
found in different bat species that are distributed in geo- 
graphically distant regions, suggesting a low degree of host 
restriction for coronavirus in those bat populations. 

In contrast to the enormous diversity of CoV genomes 
found in Old World bats [12, 25], in this study and in 
several others concerning the CoV detection in New World 
bats, only alphacoronaviruses were detected [2, 3, 17, 18]. 
Based on these results, it has been hypothesized that CoVs 
found in New World bats are less diverse than those 
detected among Old World bats [26]. 

In this initial study, samples were restricted in location 
and variety of bat species, and we found only alphacoro- 
naviruses. Such findings do not reflect data on incidence or 
prevalence of such infections in bat populations. However, 
one cannot exclude the possibility that a greater diversity 
may become apparent in Brazilian bats as long as larger 
numbers of samples from a wider spectrum of species are 
examined. 

To our knowledge, this is the first report of CoV 
detection in feces from presumably healthy insectivorous 
bats in Brazil. However, it is very likely that other bat 
species might also be infected with similar viruses. Addi- 
tional studies with larger numbers of bats and bat species, 
as well as the continued vigilance on the occurrence of 
viral infections in bats over the years is required to follow 
the evolution of bat coronaviruses in its interactions with 
the different bat host species. In addition, the detection of 
CoVs in Brazilian bat populations in close proximity to 
human inhabitants may represent a risk to human health. 
Our findings point to the need to identify the prevalence of 
CoVs in Brazilian bats, to perform risk assessment studies 
and continued surveillance of coronavirus infections from 
both urban and rural environments. 
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